Yates analysis

Full- and fractional-factorial designs are common in designed experiments for engineering and scientific applications. In these designs, each factor is assigned two levels. These are typically called the low and high levels. For computational purposes, the factors are scaled so that the low level is assigned a value of -1 and the high level is assigned a value of +1. These are also commonly referred to as "-" and "+".

A full factorial design contains all possible combinations of low/high levels for all the factors. A fractional factorial design contains a carefully chosen subset of these combinations. The criterion for choosing the subsets is discussed in detail in the fractional factorial designs article.

Formalized by Frank Yates, a Yates analysis exploits the special structure of these designs to generate least squares estimates for factor effects for all factors and all relevant interactions. The Yates analysis can be used to answer the following questions:

What is the ranked list of factors?
What is the goodness-of-fit (as measured by the residual standard deviation) for the various models?

The mathematical details of the Yates analysis are given in chapter 10 of Box, Hunter, and Hunter (1978).

The Yates analysis is typically complemented by a number of graphical techniques such as the dex mean plot and the dex contour plot ("dex" stands for "design of experiments").

1 Yates order
2 Output
3 Parameter estimates as terms are added
4 Model selection and validation
5 Graphical presentation
6 Related Techniques
7 External links
8 References

Yates order

Before performing a Yates analysis, the data should be arranged in "Yates order". That is, given k factors, the k^th column consists of 2^{(k - 1)} minus signs (i.e., the low level of the factor) followed by 2^{(k - 1)} plus signs (i.e., the high level of the factor). For example, for a full factorial design with three factors, the design matrix is

$- - -$

$%2B - -$

$- %2B -$

$%2B %2B -$

$- - %2B$

$%2B - %2B$

$- %2B %2B$

$%2B %2B %2B$

Determining the Yates order for fractional factorial designs requires knowledge of the confounding structure of the fractional factorial design.

Output

A Yates analysis generates the following output.

A factor identifier (from Yates order). The specific identifier will vary depending on the program used to generate the Yates analysis. Dataplot, for example, uses the following for a 3-factor model.

1 = factor 1

2 = factor 2

3 = factor 3

12 = interaction of factor 1 and factor 2

13 = interaction of factor 1 and factor 3

23 = interaction of factor 2 and factor 3

123 = interaction of factors 1, 2, and 3

A ranked list of important factors. That is, least squares estimated factor effects ordered from largest in magnitude (most significant) to smallest in magnitude (least significant).
A t-value for the individual factor effect estimates. The t-value is computed as

$t = \frac{e}{s_e}$

where e is the estimated factor effect and s_e is the standard deviation of the estimated factor effect.

The residual standard deviation that results from the model with the single term only. That is, the residual standard deviation from the model

$\textrm{response} = \textrm{constant} %2B 0.5 X_i$

where X_i is the estimate of the i^th factor or interaction effect.

The cumulative residual standard deviation that results from the model using the current term plus all terms preceding that term. That is,

$\textrm{response} = \textrm{constant} %2B 0.5 \mathrm{(all\ effect\ estimates\ down\ to\ and\ including\ the\ effect\ of\ interest)}$

This consists of a monotonically decreasing set of residual standard deviations (indicating a better fit as the number of terms in the model increases). The first cumulative residual standard deviation is for the model

$\textrm{response} = \textrm{constant}$

where the constant is the overall mean of the response variable. The last cumulative residual standard deviation is for the model

$\textrm{response} = \textrm{constant} %2B 0.5 \mathrm{(all\ factor\ and\ interaction\ estimates)}$

This last model will have a residual standard deviation of zero.

Parameter estimates as terms are added

In most cases of least squares fitting, the model coefficients for previously added terms change depending on what was successively added. For example, the X₁ coefficient might change depending on whether or not an X₂ term was included in the model. This is not the case when the design is orthogonal, as is a 2³ full factorial design. For orthogonal designs, the estimates for the previously included terms do not change as additional terms are added. This means the ranked list of effect estimates simultaneously serves as the least squares coefficient estimates for progressively more complicated models.

Model selection and validation

From the above Yates output, one can define the potential models from the Yates analysis. An important component of a Yates analysis is selecting the best model from the available potential models. The above step lists all the potential models. From this list, we want to select the most appropriate model. This requires balancing the following two goals.

We want the model to include all important factors.
We want the model to be parsimonious. That is, the model should be as simple as possible.

In short, we want our model to include all the important factors and interactions and to omit the unimportant factors and interactions. Note that the residual standard deviation alone is insufficient for determining the most appropriate model as it will always be decreased by adding additional factors. Instead, seven criteria are utilized to define important factors. These seven criteria are not all equally important, nor will they yield identical subsets, in which case a consensus subset or a weighted consensus subset must be extracted. In practice, some of these criteria may not apply in all situations, and some analysts may have additional criteria. These criteria are given as useful guidelines. Mosts analysts will focus on those criteria that they find most useful.

Practical significance of effects
Order of magnitude of effects
Statistical significance of effects
Probability plots of effects
Youden plot of averages
Practical significance of residual standard deviation
Statistical significance of residual standard deviation

The first four criteria focus on effect sizes with three numeric criteria and one graphical criterion. The fifth criterion focuses on averages. The last two criteria focus on the residual standard deviation of the model. Once a tentative model has been selected, the error term should follow the assumptions for a univariate measurement process. That is, the model should be validated by analyzing the residuals.

Graphical presentation

Some analysts may prefer a more graphical presentation of the Yates results. In particular, the following plots may be useful:

Ordered data plot
Ordered absolute effects plot
Cumulative residual standard deviation plot

Related Techniques

Multi-factor analysis of variance
Dex mean plot
Block plot
Dex contour plot

External links

Yates Analysis

References

Box, G. E. P.; Hunter, W. G., and Hunter, J. S. (1978). Statistics for Experimenters: An Introduction to Design, Data Analysis, and Model Building. John Wiley and Sons. ISBN 0471093157.

This article incorporates public domain material from websites or documents of the National Institute of Standards and Technology.